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ABSTRACT 

The distributions and abundances of small organics in protoplanetary disks are potentially powerful 
probes of disk physics and chemistry. HNC is a common probe of dense interstellar regions and 
the target of this study. We use the Submillimeter Array (SMA) to observe HNC 3-2 towards the 
protoplanetary disks around the T Tauri star TW Hya and the Herbig Ae star HD 163296. HNC 
is detected toward both disks, constituting the first spatially resolved observations of HNC in disks. 

We also present SMA observations of HCN 3-2, and IRAM 30m observations of HCN and HNC 1-0 
toward HD 163296. The disk-averaged HNC/HCN emission ratio is 0.1-0.2 toward both disks. Toward 
TW Hya, the HNC emission is confined to a ring. The varying HNC abundance in the TW Hya disk 
demonstrates that HNC chemistry is strongly linked to the disk physical structure. In particular, the 
inner rim of the HNC ring can be explained by efficient destruction of HNC at elevated temperatures, 
similar to what is observed in the ISM. To realize the full potential of HNC as a disk tracer requires, 
however, a combination of high SNR spatially resolved observations of HNC and HCN, and disk 
specific HNC chemical modeling. 
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1. INTRODUCTION 

The disks around young stars contain the gas 
and dust from which planets form . From theory 
( Aikawa et al.lll996l[B ergin et ^l2003fl and observations 
( Rosenfeld et al.ll2012t iQi et al.ll2013all . disks are char¬ 
acterized by steep temperature and radiation gradients 
both radially and v ertically. These disk s t ructures affect 
plane t formation (jTakenchi et al.l 119961: iKretke fc LinI 
l2012t) . motivating the quest for new emission-line di¬ 
agnostics that can trace physical structure gradients in 
disks. 

This study focuses on HNC, a common high density 
gas tracer in the ISM. At low temperatures, HNC forms 
efficiently together with HCN through ion-molecule re¬ 
actions. HNC and HCN are observed in a variety 
of different sources, but their relative abundances can 
vary by more than order of magnitude, which has 
been related to the physical environment where the 
emission originates. In pre-stellar and protostellar en¬ 
vironments, the observed HNC/HCN abundance ra¬ 
tio is ~1 at 10 K and decreas e s as temperature in ¬ 
creases fe.g. ISchilke et al.l (1199^: iPratau e t al.l (I19 97D: 
ILiszt fc Lucasl (l2OM U J0rgensenl 1)200411 : iPadovani et al l 
( 201111 : iJin et al.l (|2015fl l. suggesting that thermal reac¬ 
tion barriers pl ay an important ro le in controlling this 
ratio. Work bv iLoison et al.l (|2014ll has shown that the 
HNC/HCN ratio also depends on the elemental abun¬ 
dances of carbon at low temperatures, due to destruc¬ 
tive reactions with atomic carbon. In less shielded re¬ 
gions, th e HNC/HCN ratio als o depends on the radiation 
present (|Meiierink et al.l [200711 . 


dgraninger@cfa.harvard.edu 


Nearby, protoplanetary disks, such as the one sur¬ 
rounding TW Hya, provide an interesting test bed for 
evaluating the relative impact of different environmen¬ 
tal parameters on the HNC chemistry, due to their rela¬ 
tively well characterized geometries and density, temper¬ 
ature and radiation structures. The HNC distribution 
in disks is interesting in its own right. HNC and HCN 
are two of ver y few nitrogen-bearing molecules detected 
toward disks (IPutrey et al.l 119971 l2007t iChapillon et al.l 
120121 : lOberg et al.l 1201511 and their distributions may 
therefore proyide some of the best ayailable constraints 
on the nitrogen reseryoir. Cyanides and iso-cyanides 
are also important for s ome prebiotic chemistry scenar¬ 
ios (|Powner et al.l [200^ . Obserying HNC in disks has 
proyen challenging, howeyer. The intrinsic disk scales 
are small and the molecular line fluxes are weak com¬ 
pared to interstellar clouds and protostars. To date, 
HNC has only been detected in a single disk, DM Tau, 
(iDutrey et al.l 1199'^ and no spatially resolyed obserya- 
tions exist. 

In this letter, we report the first interferometric de¬ 
tection of HNC in tw o disks, a r ound the T Tauri star 
TW Hya (D = 54 pc IQi et al.l (1200^1 and Herbig Ae 
star HD 163296 (D = 122 nc. lQi 0470 (120011 1. using the 
SMA. We combine these obseryations with new and pre- 
yious SMA obseryations of HCN toward TW Hya and 
HD 163296, and new IRAM 30m obseryations of HNC 
and HCN 1-0 toward HD 163296, with the aim to con¬ 
strain the HNC/HCN chemistry in disks. 

2. OBSERVATIONS OF HNC AND HCN 

The HNC 3-2 line at 271.981 GHz was obseryed to¬ 
ward the TW Hya disk (Aa = ll'‘01™51.875", = 
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-34°42'17.155"; J2000.0) on 2013 April 5 and 2014 April 
8 with the SMA eight-antenna interferometer located on 
Manna Kea, HI. Observations were taken using the com¬ 
pact (COM) and extended (EXT) arrays configurations 
(projected baselines of 6.3-164.8 m). The gain calibrator 
was J1037-295. Calisto and Titan were used as the flux 
calibrators and the bandpass response was calibrated us¬ 
ing 3C279. The derived fluxes of J1037-295 were 0.79 Jy 
on 2013 April 5 and 0.70 Jy on 2014 April 8. 

The same HNC line was observed toward the disk 
around HD 163296 (Aa = 17^56'"21.28®, AJ = - 
21°57'22.1"; J2000.0) on 2013 April 5 with the SMA in 
the compact (COM) array configuration with projected 
baselines of 7.7-62.1 m. J1733-130 was used as the gain 
calibrator. The bandpass response was calibrated using 
BLLAC. Flux calibration was done using observations of 
Neptune. The derived fluxes of J1733-130 was 1.34 Jy. 

The HCN 3-2 line at 265.886 GHz was observed toward 
HD 163296 on 2008 April 20 and May 12 with the SMA 
in the compact-north (COM-N) array configuration with 
projected baselines of 7.8-114.3 meter. J1733-130 was 
used as the gain calibrator. The bandpass response was 
calibrated using observations of 3C273. Flux calibration 
was done using observations of Titan and Callisto. The 
derived fluxes of J1733-130 were 2.50 Jy on 2008 April 
20 and 2.13 Jy on 2008 May 12. 

All data were phase- and amplitude-calibrated using 
the MIR software packagcQ. Continuum and spectral line 
maps were generated and CLEANed using the MIRIAD 
software package. Images were constructed from the TW 
Hya data using a Briggs (robust) weighting factor of 0 
and in HD 163296 data using 2.0 and 0.3 for HCN and 
HNC, respectively. The resulting synthesized beam for 
HNC in TW Hya is 1"8 x T.'O. For HD 163296, the syn¬ 
thesized beams for HCN and HNC are 370 x 179 and 
475 X 275, respectively. 

HCN 1-0 and HNC 1-0 lines were observed as a part 
of a single-dish 3 mm line survey of the disk around HD 
163296 on 2012 September 13 - 16 with the IRAM 30m 
telescope on Pico Veleta using the EMIR 90GHz receiver 
and the Fourier transform spectrometer (FTS) backend 
in average summer weather (pwv 4-8 mm). Two spec¬ 
tral settings were used with the four sidebands covering 
84.0 to 115.8 GHz with a resolution of 200 kHz (^0.5 
- 0.6 km/s). Pointing was checked every 1-2 hrs and 
focusing every 4 hrs. The data was acquired in the wob¬ 
bler switching modes, with typical throws of 60-120” and 
an on-source integration time of ^12-15 hrs per spectral 
setting. Spectra were baseline subtracted and averaged 
using CLAS^. To convert from antenna temperature, 
T*, to main beam temperature, T^b, forward and beam 
efficiencies of 0.95 and 0.81 were applied. 

3. RESULTS 
3.1. Observations 

HNC 3-2 emission is clearly detected and resolved to¬ 
ward TW Hya. Figure [1] displays the 1.1 mm continuum 
and HNC 3-2 SMA observations, combining the com¬ 
pact and extended data. The resolved integrated con- 
tinu um flux is 0.78 Jy which agrees with previous stud¬ 
ies ()Qi et al.ll2013all . The velocity-integrated HNC flux 

^ http://www.cfa.harvard.edu/~cqi/mircook.html 

^ http://www.iram.fr/IRAMFR/GILDAS 


peaks at 0.53 Jy km s“^ beam“^. The integrated HNC 
emission peaks are offset to the east and west of the con¬ 
tinuum phase center, suggestive of a ring-like emission 
structure within the disk. In the HNC image the west¬ 
ern flux peak is almost twice as strong compared to the 
eastern one, corresponding to ~ Scr. Asymmetries of this 
magnitude may be artifacts of CLEANing process and 
based on visual inspection alone we are agnostic about 
its significance. We evaluate this aspect of the the data 
when modeling the complex visibilities in §3.2. 
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Figure 1. The 1.1 mm continuum and HNC 3-2 emission (up¬ 
per panels) from TW Hya using the compact -|- extended config¬ 
urations of the SMA. Continuum contours are [10,20,40,...] mjy 
and the contours for HNC are [4,6,8 ](t with Icr = 0.064 Jy km 
beam~^. The lower left panel displays the spectrum obtained from 
the HNC 3-2 image. The lower right panel displays the deprojected, 
azimuthally averaged flux as a function of radius. 

To isolate the radial flux variation in the HNC emis¬ 
sion, Fig. [1] (lower right panel) shows the deporjected 
azimuthally averaged HNC line flux, usin g an inclination 
angle of 6° and positio n angle of 155° (iQi et al.l 12004 
I Williams fc: Cieza 1201 111 . The resulting emission profile 
reaches a maximum at ~0.5”, indicating the presence of 
a small HNC hole at the center of the TW Hya disk. 
This agrees with the separation between the two 
flux peaks in the HNC integrated emission image. The 
depth and size of the hole may both be affected by the 
image asymmetry, however. In §3.2 we evaluate the ev¬ 
idence for an HNC ring and quantify its location based 
on an analysis in the (it, ii)-plane. 

The HNC spectrum in Fig. [T] was extracted using a 
275 X 275 box. The integrated line flux is 1.2 Jy km s“^ 
beam“^ (~12 (t). Using the the inte grated HNC flux from 
this study and the HCN flux from Qi et al.l (|2008D . the 
disk-averaged HNC/HCN emission ratio is 0.18 ± 0.06, 
where the error is derived using the Icr rms for the two 
observations. 

HNC is also detected toward HD 163296, but with a 
lower SNR. Figure 2 displays the 1.1 mm continuum, 
HNC, and HCN emission toward HD 163296 as well as 
the spectra of the molecules and the velocity structure of 
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Figure 2. 1.1 mm continuum (top left, contours=[20,40,80,...] mjy), integrated HNC 3—2 (top middle) and HCN 3-2 (top right) fluxes 
toward the disk around HD 162936 together with extracted spectra (bottom left) and integrated emission in two velocity bins (bottom 
middle and right). The HNC and HCN integrated flux contours are [2,3]o- (HNC) and [3,5,7,..]a- (HCN) with Icr =0.21 and 0.14 Jy km s“^ 
beam~^, respectively. The HNC and HCN contours in the two velocity bins are [3,4 ](t (HNC) and [8,10,12,..]cr (HCN). 

HNC and HCN. The velocity integrated emission maps 
yields a maximum HNC and HCN flux of 0.71 Jy km s“^ 
beam“^ and 2.71 Jy km s“^ beam~^, respectively. The 
HNC emission toward the HD 163296 disk is suggestive 
of a ring viewed at relatively high inclination, but the 
low SNR inhibits any further spatial analysis. 

Spectra were extracted using an elliptical mask, which 
increased the SNR compared to a square mask. The 
total integrated flux from these spectra is 0.66 Jy km/s 
beam“^ (~4 (t) and 4.47 Jy km/s beam“^ (^9.5 (t) for 
HNC and HCN, respectively. This yields an HNC/HCN 
ratio of 0.14 ± 0.05, similar to TW Hya. 

The IRAM 30m observations of HNC and HCN 1-0 
toward HD 163296 should be more sensitive to the coldest 
disk material (Fig. [3]). The integrated intensities of HNC 
and HCN were found to be <0.033 K km s“^ (3tT upper 
limit) and 0.104 K km s“^, respectively. This yields an 
upper limit on the HNC/HCN ratio of <0.3, consistent 
with the SMA data. 

3.2. HNC line modeling toward TW Hya 

To constrain the spatial extent of the HNC e mission in 
TW Hya, we follow the procedure outlined in IQi et ^ 

(|2008D : lOberg et al.l (|2012l I2015D to extract molecular 
column density profiles. First, parametric HNC abun¬ 
dance models are defined with respect to the TW Hya 
disk te mperature and density structures from IQi et all 
(|2013bD . The line emission is modeled using the two- 
dimensional no n-LTE Monte Carlo radiative tr ansfer 
code RATRAN (|Hogerheiide fc van der Ta3l2000D using 



Figure 3. IRAM 30m spectra of HNC and HCN 1-0 in HD 
163296. The dashed line at 5.8 km/s is the systemic velocity of 
HD 163296. 

collisional cross sections from iDumouchel et aj ] (IMni) 
The H NC abundance models follow the form in IQi et all 
(|2008D , which assumes a molecular layer with a vertically 
constant abundance, bound by surface (ds) and midplane 
(cm) boundaries defined with respect to the vertically 
integrated hydrogen column in the unit of 1.59x10^^ 
cm~^. We set t he HNC boundaries using the values from 
IQi et al.l I12008D for HCN (logcrs=-0.5, log(Tm=1.5). The 
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Figure 4. Comparison of HNC 3-2 data (upper left panel) and models with an outer (Ro) and hole (R^,) radius. The contours are [3,4,5]cr, 
based on the observed RMS noise. Images are generated using the EXT data only to display the small-scale emission structures. 


radial distribution of HNC was modeled using an abun¬ 
dance power-law with an outer radius (Ro), hole radius 
(R/i), and power-law index of -1.5. The parameter Ro 
was varied between 35 and 60 AU and R^ was varied 
between 20 and 50 AU. The minimum value was cal¬ 
culated from the modeled and observed HNC emission 
in the (u, r!)-plane using the data from the compact and 
extended configuration. 

Figure |4] shows examples of emission models and ob¬ 
servations in the extended configuration. The best fit 
model based on the value is a 10 AU wide ring with 
a Ro of 45 AU and R^ of 35 AU (Figure 0, shown in 
the lower left panel of Figure [H This model results in 
no significant residuals when subtracted from the data 
in the {u, n)-plane and then imaged. The observed HNC 
azimuthal asymmetry toward TW Hya is thus not sig- 
nihcant. The presence of a HNC ring in the TW Hya 
disk is by contrast a statistically signihcant result. Fig¬ 
ure [5] shows the models that fit the data within 2 and 3 (t. 
Within 2cr, the hole size is 28-45 AU (3cr: 26-47 AU) and 


the outer edge is at 37-55 AU (3cr: <57 AU). 

4. DISCUSSION 

In this study, we detect and spatially resolve HNC 
emission toward TW Hya and HD 163296. From these 
observations we extracted disk averaged HNC/HCN 
emission ratios, and the inner and outer HNC emission 
radii in the TW Hya disk. Here, we will discuss the 
constraints that these observations provide on the HNC 
chemistry and the disk structure. 

4.1. HNC/HCN Emission Ratios 

The disk-averaged HNC/HCN emission ratio is 0.1 - 
0.2 toward both T W Hya and HD 163296. Based on 
ISM observa tions (ISch ilke et ^ 119921) and astrochem- 
istry models (|Graninger et al.ll2014fl . the HNC/HCN ra¬ 
tio can be used to constrain the characteristic tempera¬ 
ture of the emission region because of efficient destruc¬ 
tion of HNC by hydrogen atoms at temperatures >25 K. 
However, to apply this result straightforwardly requires 
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Figure 5. contour plot for the best fit models. The contours 

are 2 and 3cr. The 3(7 contour is in red and the red cross marks 
the best-fit model. The black crosses mark the values of and 
Ro modeled and the shaded region indicates where R;i>Ro. 


that HNC and HCN are co-spatial, that there is a sim¬ 
ple relationship between emission and abundance ratios 
of HNC and HCN, and that temperature dominates the 
chemistry, none of which clearly hold in the case of these 
disks. 

First, HNC and HCN are not co-spatial - some 
HCN is present toward the center of the disk, where 
HNC is absent. This may be partially mitigated by 
averaging the HNC and HCN emission only across 
disk areas where HCN and HNC emission are respec¬ 
tively detected. Second, the relationship between line 
emission and abundance ratios is complicated by high 
HCN line opacitie s and different HCN and HNC exci¬ 
tation properties (jSarrasin et al.l [2?)T(TI1 . We used the 
large velocity gradient ra diative transfer code, RADEX 
(jvan der Tak et ;i I]|2007D , wi t h the collisional rate con- 
stants from iDumouchel et al.l (1201011 to evaluate relation¬ 
ship between the HNC/HCN line emission and abun¬ 
dance ratios using temperat ures and densities c haracter- 
istic of HCN disk emissionjWalshe^jb 12010( 1 and the 
HCN column density from iQi et al.l (|200S 1. We find that 
the line and abundance ratios can differ by up to an order 
of magnitude, demonstrating that more detailed model¬ 
ing is required to extract HNC/HCN abundance ratios 
from the observed line emission. 

Even if abundance ratios are successfully extracted, it 
may be difficult to infer a characteristic emission region 
temperature. Other environmental characteristics than 
temperature may regulate th e HNC/HCN abundanc e ra¬ 
tio in disks. The model by iGraninger et al.l (1201411 ne¬ 
glects the influence of radiation, which is reasonable in 
cloud cores, but may not b e a good appro x imatio n in 
disks (|Meiierink et al.l [2007( 1. iLoison et al.l (j2014l l has 
also shown that the HNC/HCN ratio depends on the 
C/0 ratio, which may be elevated in disks. In the ab¬ 
sence of detailed disk chemistry models that incorporate 
the complete HNC/HCN chemistry, disk-averaged line 
ratios are thus not useful probes of disk properties. 


4.2. The HNC inner and outer radii in the TW Hya 

disk 

In the disk toward TW Hya, there is a lack of HNC 
emission interior to 28-45 AU, co rresponding to te mper- 
atures above 21-13 K using the lOi et al.l (|2013bll tem¬ 
perature pro^e. HCN does not have an emission hole 
ijOi et ahlboOSH . The observed absence of HNC close to 
the star is c onsistent with tempera ture dependent HNC 
destruction (jCraninger et al.ll2014r i. It is difficult to ex¬ 
plain the hole with the increased UV flux present at 
smaller radii, since HNC and HCN are t hought to have 
similar photodissociation cros s-sections (|Roberge et al.l 
(19911 : l^n Dishoeck et al.ll200fil l. 

HNC also has an outer emission radius between 35-55 
AU, with the best fit at 45 AU, substantially smaller than 
the observe d radii of othe r gases, including CO, HCO'*' 
and HCN (|Qi et al.l l2008fl . This outer cutoff in HNC 
emission was unexpected. It coincides with the onset 
of C 2 H emission, wh ich is proposed to b e a UV pho¬ 
tochemistry product (jKastner et al.l[2?)T^ . and is close 
to the radius where t he dust density sharply decreases 
(j Andrews et al.l [2012( 1. These spatial coincidences sug¬ 
gest three potential explanations for the outer cutoff of 
HNC emission: 1. UV photodissociation of HNC due to 
deeper penetration of stellar and interstellar UV photons 
in the low-column density outer disk, 2. an increased UV 
photodissociation of H 2 and carbon-bearing molecules 
accelerating the HNC+H and H NC+C destruction pro¬ 
cesses at a given temperature (|Graninger et al.l (2014 
ILoison et al.ll20l4l . and 3. grain-gas thermal decoupling 
due to lower gas and grain densities, causing a temper¬ 
ature reversal in the outer disk and thus the increase 
in thermally controlled HNC destruction (i.e. a similar 
process that sets the inner edge). 

The first two explanations both depend on a more 
efficient UV photochemistry exterior to 45 AU, consis¬ 
tent with the onset of C ^H emission at the same ra¬ 
dius (jKastner et al.ll2?)T^ . and signs of efficient UV pho¬ 
tochemistry in the outer low-density regions of other 
disks (lOberg et al.l 120111) . Flared disks can intercept 
stellar UV radiation out to larger radii, and Lya emis¬ 
sion is efficien tly down-scattere d through disk towards 
the midplane (jFogel et al.ll20lj] ). Additionally, the in¬ 
terstellar radiation field will contribute a non-negligible 
amou nt of UV radiation to the outer disk (| Her gin et al.l 
l2003f) . Efficient UV penetration at large disk radii may 
increase the HNC photodissociation rate, reducing its 
abundance. This scenario could be tested by compar¬ 
ing HNC and HCN emission profiles. HNC and HCN are 
believed to have simi l ar photodissociation cross -sections 
(jRoberge et al.l [T99H Ivan Dishoeck et all 1200611 . and if 
their outer emission radii are UV controlled, they should 
be the sa me. The best-fit outer HCN emission radius is 
'^lOO AU (|Qi et al.l 120081) . twice that of HNC, but this 
value is not well enough constrained to exclude direct 
photodissociation of HNC as an explanation for the HNC 
outer radius. 

The second and third proposed explanations for the 
HNC outer cut-off should produce a smaller HNC than 
HCN emission radius, due to an increased HNC-to-HCN 
conversion rate in the outer disk. In the second scenario 
and increased UV penetration increases the atomic H 
and C abundance through H 2 and CO photodissociation 
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(|Akimkin et alJl2013D. speeding up HNC iso merization at 
a given temperature ([Graninger et aljr2014f) . In the third 
scenario, a lower density in the outer disk reduces gas- 
grain thermal coupling sufficiently that direct gas heat¬ 
ing by the interstellar radiation field becomes important. 
This heating could result in a radial temperature gradi¬ 
ent reversal outside of 45 AU, again increasing the HNC 
isomerization rate. 

Higher spatial resolution observations with ALMA 
should be able to distinguish between the three scenar¬ 
ios. The relative sizes of the HNC and HCN outer radii 
should differentiate between ^1 on one hand and #2 and 
^3 on the other. If ^1 is ruled out, the radial temper¬ 
ature profile of HCN beyond 45 AU should be able to 
differentiate between ^2 and #3. Whichever scenario is 
confirmed, spatially resolved HNC/HCN emission mod¬ 
eled using non-LTE techniques clearly has a large poten¬ 
tial to probe interesting disk properties. 
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